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Abstract. We study the fundamental problem of two gas species whose molecules collide as hard spheres 
in the presence of a flat boundary and with dependence on only one space dimension. More specifically 
the steady linear problem considered is the one arising when the second gas dominates as a flow moving 
towards the boundary with constant microscopic velocity (and hence zero temperature). The boundary 
condition adopted consists of prescribing the outgoing velocity distribution of the first gas at the boundary. 
It is discovered that the presence of the boundary under general assumptions on the outgoing distribu- 
tion ensures the convergence of a series of path integrals resulting in a convenient representation for the 
distribution of the velocities of the molecules returning at the boundary. 



We begin by noting that the Boltzmann equation with boundary conditions has been previously con- 
sidered in a number of works, for example the reader may consult the work by Bardos, Caflisch, and 
Nicolaenko PQ, the zero-temperature paper by Caflisch [3] where behavior in the interior is studied in a 
perturbative setting with regard to applications to strong shock waves, as well as the works [8] and [12] 
among many others. Our approach differs substantially from the previous works since we investigate the 
structure of the solution in terms of contributions from the individual particle paths and we look for an 
exact description of the solution at the boundary in the velocity variable. Excellent reviews of kinetic 
theory in general are [6] and [Uj as well as [T5] . 

In our problem we consider molecules emitted from a flat boundary which interact as hard spheres with 
the molecules of a second gas which have zero temperature and penetrate into the boundary. We will be 
interested in the distribution of the velocities of the emitted molecules at the time of their return to the 
boundary. The interactions between the emitted molecules themselves are ignored and so a linear problem 
is studied. Accounting for the nonlinear interactions would be an important continuation of this work. We 
must also point out that the zero temperature distribution of the background is of course an idealization. 
However it is reasonable to expect that many of the features studied here will be still present in the case 
of the more natural Gaussian distribution of the background (see the concluding remarks). 

The equation at hand is time independent, valid for x > 0, and takes the form 



where /(£, x) with £ = £ 2 , £ 3 ) is the unknown velocity distribution and p5(£ — c) with c = (— c, 0, 0) for 
c > is the background distribution where p is the number density of the background. So the background 
molecules enter the boundary and their distribution is homogeneous in the interior. Here Q is the Boltz- 
mann collision operator. To simplify the presentation let us assume that both particle species have the 
same molecular mass, normalized to 1. The case when the masses differ can be treated similarly. We will 
focus on determining the value = /(£,0) for f 1 < 0. After writing out Q as the difference between 

the gain and loss terms due to collisions we have 



1. Introduction and Statement of the Result 



ed x f(t,x) = Q(f,pS(i;-c)) 

/tt,0)=/ + (0far€ 1 >0 





/(£ - ((£ - £*) • n)n)5(& + ((£ - £.) • n)n - c)|(£ - £*) • n\dnd£* 
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Figure 1 . The plane of influence 
= pa 2 Kf - P a 2 K\i - c\f. 

where n is the unit vector at which the collision occurs and S + is the half sphere (£ — £*) • n > 0, a is 
the sphere radius. For our purpose we will take pa 2 = 1. This is actually no restriction as explained in 
the remark on scaling after the statement of the theorem. We will give a more convenient expression for 
Kf(£o), but first to elucidate it we will determine the region containing particles £i which can influence 
the region near £ after experiencing a collision at the direction of n. In other words given £ = (£o > Co > Co) 
we need to determine the set {n, £i = (CijCi>Ci)} such that 

Co = Ci - ((6 - c) • n)n. 

In the coordinate system n,n let P„ and P^ denote the corresponding projections. Since the effect of 
a collision is to exchange the momentum along n, the new momentum along n of the particle hit by a 
background particle is always P n (c), i.e. we see that P n (c) = P n (Co)- This is satisfied iff n _L £o ~ c which 
determines n. But then we must also have -P^(Ci) = -PrMCo)- Therefore the particles with velocity £i that 
can influence the velocity £o lie along the plane L c ^ passing through £ and perpendicular to £ — c (see Fig. 
1). In fact it can be easily shown that the interior of the sphere C^ on the figure is the region influenced 
by £o after an arbitrary number of collisions and the exterior is the region that influences £o- This answers 
our question. So clearly K will be given by an integral along this plane but we need to determine the 
correct weight. Although a direct calculation is of course possible, for brevity to obtain the exact formula 
we will make use of Torsten Carleman's well-known representation ([4], [5]) for the gain term, namely: 

JR 3 JS+ 

JR3 I? - SI I JL ilX 

where L^^ x denotes the plane perpendicular to £ — £i and passing through £. For completeness the proof of 
this representation is briefly indicated in the Appendix. Specializing the above to the case which interests 
us, namely <?(£) = <$(£ — c) we obtain that 

Kffa) = I /(CiRi- 

l?o - c| J Lc io 

In what follows we will omit c from the subscript of the plane since c is fixed. 
We are going to prove the following 

Theorem 1.1. Let there be given an outgoing (from the boundary) distribution / + (£) satisfying < / + < 
M which is supported in the compact region t; 1 > and |£ — c| < R. Then for the incoming (to the 
boundary) distribution denoted by / _ (£o) an d defined for £q < we have the following: 
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i) An explicit series with positive terms with factorial convergence exists such that 

oo 

r(£o) = ]T/n(£o) 

71=1 

where the n-th term represents the contribution of particles returning after n collisions, and specifically the 
rate of convergence is such that if lZ n denotes the remainder term then the mass flux due to it satisfies 

MR 5 
—ne 



So' 



whenever it is satisfied that 



n > 30(log— + 10). 

c 



We see that the number of terms describing the solution grows linearly in log — . 

ii) If we use the notation ro = |£o — c | the solution /^~(£o) is identically for r® > R and it possesses a 
singularity at ro = 0. The specific nature of the singularity is such that for ro < c/2 and any S > we have 

MR 3 1 

were C$ is an explicit constant whose dependence on c and R is indicated in the proof. 

We remark that although the above bound already indicates fast convergence, a more careful analysis 
may improve notably the constants given. But for us it is most important that these constants are explicit 
and that in our proof we establish the mechanism leading to this convergence. The estimates indicated 
above present integrable functions and the first terms of the series mentioned above contain essentially all 
the mass flux of the solution. The mass flux is the natural measure of the distribution at the boundary 
(as is made precise for example in |7]). The results of this theorem are in contrast to other situations in 
kinetic theory where such expansions cannot be expected to converge with any reasonable rate. In these 
other situations the lack of fast convergence is due to the fact that (fluid) components of the solution cor- 
responding to the kernel elements of the collision operator (and representing thermodynamic equilibrium) 
propagate in the interior of the domain of interest (in our case K — v has only one kernel element which is 
the Dirac delta and our linear equation only has conservation of mass) . In our problem the characteristics 
do not propagate inside the domain so these fluid components are irrelevant and we are studying a purely 
kinetic phenomenon. To study fluid components propagating in the interior (at least in the case of Gauss- 
ian background where 1? space methods are applicable) understanding of the spectrum and the use of 
perturbation theory are most natural (see for instance [llj). Although we do not investigate this here, it is 
reasonable to expect that a description similar to ours will also apply to the kinetic part of the solution in 
the case of the so called linearized collision operator (when a perturbation of the background is considered 
so energy and momentum are also conserved) but in such a case depending on the speed of the background 
and the conditions at infinity a fluid component of the solution may also appear as a remainder term in 
the above series. 

The behavior ro — > presented by the above formula is essentially optimal, in fact except for the 6- 
correction it is already displayed by the first term of the series. In the limit c — > we have Cs/c — > oo. This 
reflects the fact that when the background velocity vanishes the linearized problem does not posses a finite 
steady solution with zero flux at infinity. The equilibrium in the case c sa is due to the particles emitted 
from the wall interacting with themselves which is a nonlinear effect, or alternatively a linearization is 
necessary around certain density which does not vanish away from the wall. 

We hope that the understanding gained from this problem may be useful in studying other boundary 
conditions, as for example the idealized situation of specular reflection £ 2 , £ 3 ) = /(— £ 2 , £ 3 ) at 
x = 0. In this problem it would be quite interesting to know the momentum flux / K s^ (s ) /(C)^ a * 
x = which gives the force exerted on the boundary by the background through the interaction with the 
particles trapped near the boundary. 



4 



ALEXANDER SOTIROV 



Now we briefly discuss the scaling in this problem. There is no macroscopic reference length in the x 
variable, so a Knudsen number cannot be specified. Furthermore since we are taking x = as the region 
of attention we cannot expect our solution to be influenced by the dimensional factor pa 2 . This is in fact 
exactly so: the above theorem would hold unchanged even if there was such a factor in the right hand side 
of the original equation. The reason is that we are not interested in how long or how far a particle had to 
travel before return, we only care for the velocity at return. A proof of this is noted after equation (|2.1|) . 

2. The Particle Path Representation and Proof of the Result 

To prove the above result we will consider the contributions of the different paths a particle coming 
out of the boundary can travel until it returns to the wall as an incoming particle. More specifically 
we will distinguish these paths according to the number of collisions experienced before return. An easy 
approach to this representation comes from the Laplace transform formalism. Let £/(£, z) denote the 
Laplace transform in the x variable. We have: 

eizcif) - m,o)] = c(kj) - v c(j) => 

/(£, x) = C- 1 {[ze -K + v]- x e f(Z, 0)) 



or if we switch to the more convenient Fourier notation 



2m J -too 

convenient Fourier notation 

= ^ J e ikx [ike--K+u]- i emo)dk 



Now if we write 

[ik^ 1 -K + v]- 1 = [iki 1 + v}- 1 [I - K[ik^ + u}- 1 ] _1 
and introduce the notation £) = K[ik^ + v]^ 1 we can write the expansion 

[/ - $(*, or 1 = i + o + o 2 + • • ■ 

and we obtain a series representation for the solution. It is the convergence of the series resulting from the 
above representation that we shall establish, provided the outgoing distribution satisfies the requirements 
in the theorem. The zeroth term in the series represents particles experiencing no collision at all after 
leaving the wall, so it will be ignored. The n — th term's contribution to the incoming distribution is 

1 roo ikx 

2ttJ_ 00 ikQ +tt|£o - c| 

1 f 00 1 1 

(2-1) iz ■ 57 „ t K- 



2tt ik$ + tt|Co - c| ikQ + 7r|6 - c| 

K ikti+*\z 2 -c\ ■ ■ ■ K lk e n +^ n - c \ ieM+ ^ )dk - 



and we will construct / as 



n=l 

It is easy to see that this decomposition corresponds to the number of collisions experienced by a particle, 
since by the convolution theorem for the Fourier transform 
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X{ SS a(^)x>0} I * I e 1 X{ S gn(5})x>0} I I e *n X{sgn(ei)x>0} 



where the \ are just indicator functions. This expression is a convolution of exponential distributions (up 
to the normalizing constants) where the exponent is l/l where / is the expected distance a particle will 
travel before it experiences a collision (l/l is the ratio of the collision frequency to the velocity). Our 
formula (|2.1j) differs from this by the additional information due to K which specifies what new velocity 
the particle is likely to acquire after the collision is experienced. With this it becomes clear that iterating 
n steps represents a path consisting of n collisions. 

To address the scaling question mentioned in the introduction, note that if we had included pa 1 then in 
(|2.1j) after changing k' — k/(pa 2 ), dk = pa 2 dk' this constant will disappear at x = as desired. At x > 
it of course doesn't disappear but will appear in an exponential factor of spatial decay. 

Now fix a set of values £o»£i> - •• )£n i n the above expression. Let it be the case that £f or . > and 
£,retj < where {fori} and {retj} are sets of indices which partition {0, 1, . . . , n} with £ = 1,2,...,?, and 
j = 1,2, ... , to where l+m = n + 1. The first correspond to particles moving forwards (away from the wall) 
and the second to returning particles (moving backwards) . Throughout its n-collision trip a particle could 
experience an arbitrary sequence of forward and backward movements. We will perform the dk integration 
with these fixed velocities. For this it is necessary to compute the following integral (distinguish \/—T from 
the index i): 

dk 



nUM&J+^/o* -c|)nr=iH fe l& t ,l+*l£re tj -c|) 
Introducing the notation 

_ K\€fon - C| AZretj - c| 

-, Oj — 



The above integral can be written as 



nu^n™!^ r dk 



* n+i nr=o i& - c i J-™ uU(ik + m) n7=ife - ' 

It is possible to evaluate the above integral explicitly. In fact its evaluation is the following Lemma 
Lemma 2.1. i) If for ai > 0,bj > we use the notation 

T i,m f „ „. , ^_ r ni i"<ii:v<. 



j^(ai,..., ai , bl ,...b m )= / . dk 

J -°° 1 li=l (*« + «i) 1 Ij=l ( fo J ~ lk ) 

then the following recursive relation holds 

jl.rn _ a l jX-l,m _j_ b m jl,m— 1 

ai + b m a t + b m 

ii) Suppose that we have k indices jujz, ■ ■ ■ ,jk such that bj 3 < B for s — 1 . . . k and we also have that 
ai > A for all i. Then we have 



i „ fc L l 



jl,m < —B K I 1 



2tt V A 



2+m 



Proof, i) Notice that 



i (ife + a; + 6„, — ifc) TT. , a,- LP™ , 

a ; J + 6 m j' ,m_1 = / - p J rffc = (a; + b m )J l ' m 

J -°° n i =iF + «.)ll F i( f 'r !i: ) 

_^ ji,m _ 0,1 jl-l,m i °m jl,m — l 

ai + b m ax + b m 



(i 
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as desired. 

ii) Notice that J l ' m {a±, . . . , a/, bi, . . . , b m ) is a symmetric function in each set of variables. So without loss 
of generality we can assume that b\ < B. 

The recurrence relation from i) tells us we can construct J l ' m by constructing a tree with a root at J l,m 
and each node having either one or two children. All the leaves of the tree contain either the expression 
Q2+fci or ^ e ex P ress i° n ai+b 2 J 1 ' 1 - This is because by elementary contour integration we have J°' s = 
J s, ° = for all s > 1. Also notice that by contour integration 

,1,1 = r Mi dk = 1 ^ = 1 & i < L B 

J-oo ( ik + ai)(~ik + h) 2vrai+6i 2tt 1 + h~ ~ 2tt 

We will prove the claim by induction on m, I, k. By the above it holds for I = 1, m = 1, k = 1. 

Case i: 6 m not in {i)^, . . . , 6j fc }. Then J 1 - 1 -" 1 and J 1 '™-- 1 both still contain arguments 6j with &j < B, 

so the induction hypothesis applies to them with k unchanged, i.e. 

jl,m = _*jl-l,m + _^_jl,m-l < max | jl-l,m jl,m- 1 1 

ai+b m ai + b m 

and the result follows by induction. 

Case 2: b m € {6^, . . . , 6^.}. Then J'.™- 1 has at least fc — 1 arguments 6^ with 6j < £? and J l ~ l m still has 
A: arguments 6j with 6j < B. So applying the induction hypothesis 



. , / i \ m+l—i -, / i \ m+l—l 



as desired. □ 



Remark 2.2. The above estimate is not optimal, but it will be by far sufficient for our purpose. Notice 
that the expressions J l,m are trivially all less than J 1 ' 1 . In particular the growth due to the (I + 4) term 
in the above is not optimal. ■ 

Let cm be the angle between and £i with respect to a pole centered at c and let = |c — Let 
li,<fii be the polar coordinates in the plane L^ i _ 1 centered at This notation is illustrated on Fig. 2. 

With this we can write (jl.ip in the more convenient form 

#/te-i) = / / lif{l h 4>i)dlid4>i. 

fi-1 Jo Jo 

Eventually we will pick the angles on as variables of integration. So we consider the vectors 

&(£o, olu 0i, a 2 , 4>2 ■ ■ ■ a n , <t> n ). 
From the definition we must restrict on £ [0, tt/2). We have that 

n-i , r 
and so 



cos cti cos ai cos a-i ■ ■ ■ cos a.; 

Notice that the angles <pi, <p2, ■ ■ ■ 4>i are not needed to express r^. This is important in what follows. Due 
to the assumption on the support of / + (the shaded area on Fig. 2) we require that c < r n < R which 
translates to 

(2.2) c < r n = - < R. 

COS Ot\ COS Q.2 ■ ■ ■ COS Oi n 

This restriction is clearly not optimal, since it includes an area larger than the support of / + . This however 
will only affect the final estimate by a fixed constant related to the ratio of the volumes of the parts of the 
shell c < r n < R on the two sides of = 0. In any case this does not affect the convergence structure that 
we study. 
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Figure 2. The support of the initial data and notation 



Let W be the set of ai in [0,7r/2) n which satisfy the above condition (|2.2j) . Also let U denote [0, 27r] n 
which is the domain of integration in the fa. We also have 



li = rj_i tancti and dU 



-dai 



With this notation, the definition of K from (jl.ip and using the Lemma we can write: 

i r°° r r i i 

2tt y^oo Ju J w ik^l + tt|£o - c| |£ - c| 



/n(&) 



1 



1 



+7r|Ci(Co,ai,0i) - c| |£i(fo, ai, ^2) - c 
1 1 



fe(^0) OH,faa2,4>2 



/ + (£n) 



«fc£i + 7r|£„ - c| |£„ - c 
<i/ n cM„_i dZ 2 (i/i 

1 



)\&(€o,ai,fa ■ ■ ■ ,a-n,fax)\ 



1 sinai 1 sina2 

r n 



da\dfa . . . d4> n dk = 
1 sin a„ 



2ir n+2 Ju Jw V r o cosai r\ cosa2 
1 1 1 



r„_i cosav, 



2 Tl 2 

cos a\ cos z a2 



1 



■ r n -i- 

■ r n ...r n 

1 > ' " da 2 da 1 



J(ai, . ..,aiM, ■■ ■ ,b m )\£n\f + (£n)da 

1 r f nr=i sm ^ Tfn 



2tt"+ 2 



, ai,bi, . . . , b m )^-f + (^ n )da n . . . da 2 da x dfa 

Tn 



w j w (nr=i c ° sQ! i) 3 ' 

where J(a 1; . . . , a;, b±, . . . , b m ) denotes the integral from Lemma 2.1. The numbers I and m here are a 
function of the specific choice £o)£i> • • -£n- Now using (|2.2p we arrive at 



(2.3) /"(&) < 



1 



7 



■>-, + :< / / r 3 I II Sina J(fl1 ' 



,U — / + (^){da i # i }? =1 . 



We will often refer to this representation in what follows. When we attempt to bound the above 
expression there will be a factor of (27r)™ from the integration in the fa variables and there will be also a 
factor coming from an appropriate estimate of the J-contribution. We want to control these factors. To 
obtain a good estimate we will split the set W in two parts W = W\ + W^. On the set W\ the J-term 
inside the integral will posses geometric decay with n, and on the other hand after a change of variables 



s 
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the measure of the set of integration corresponding to W2 will decay with n sufficiently fast even when r 
is small. We need to do this in order to control /~ in such a way that their sum remains an integrable 
function; if we did not require this, a simpler estimate would still show factorial convergence, but only 
pointwise in r . Our extra effort is justified since /~ is a probability distribution and it is most natural 
in addition to pointwise estimates to establish the control of the decay of f~ in the weighted L 1 ^) norm 
/{ 1 <o Is representing the mass flux through the boundary of returning particles. We denote from 

now on this norm by L^. 

Pick < (3 < 1 and define k = [pn] + 1 where the brackets indicate the integer part. Also let D > 1 be 
a constant dependent on (3 but otherwise fixed that will be specified later. Let W\ and W2 be defined as 
follows: 

W\ : (ai, a 2 , . . . a n ) G W such that r k < c/D, 
W 2 = W — W\ i.e. W2 ■ («i, a 2 , . . . a n ) such that r k > c/D. 
Notice our definition makes sense since r\ , . . . , r k are completely determined by ct\ , . . . , a k and r$ so no 
condition on the (pi variables needs to be imposed. Now with (12. 3p in mind we write 

where in /" the region of integration in (|2.3|) with respect to the atj is W\, and that for I£ is W2. To 
evaluate the integration over these two regions we employ the following lemma. It does not focus on the 
contribution of the J-part of the integrand which will be bounded at the next stage. 

Lemma 2.3. The following evaluations hold (assuming rg < c/2): 
i) 

/ I I sin ctidctn ■ ■ ■ dai < If sin cnda n . . . da\ < — — ( log — 

Jw 1 fj{ Jw fj{ n\ c \ r Q 



ii) 



[[sma i da n ...da 1 <--^[ k _ j )(log— ) (log- 



where we recall that k = \j3n] + 1. 

Proof, i) Pass to the variables of integration = cos a, . We obtain 



/ 1 | sin ctida n . . . dai — / dx n . . . dx\ 

Jw i=1 Jn 



where we have Xi £ (0, 1] and f2 with the restriction from (|2 . 2[) is described as: 

iZ : — > X!X 2 ■ ■ ■ x n > — . 
c R 

So we need to compute Vol(fl). If we further pass to the variables yi = — logXi we get that the corre- 
sponding yi region $(0) is described as 

c " R 
: < log — < > yi < log — and yi S [0, 00). 

Furthermore since Xi — e~ Vi the Jacobean of the transformation is 

/ e- yi ... \ 

D(X1,X2, ■ ■ ■ ,X n ) 



Det 



e' V2 ... 



D(yi,V2, ■■■y n ) 

\ 

On our region of integration this Jacobean can be bounded as 

e ^1=1 < e = — 



-!),< 
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If A™ denotes the n-dimensional simplex with side x we have that 

x n 



Vol{A n x ) 

so since is the difference of two such simplices we conclude that 

VoW) < r -±[Vol(Al sJL ) - Vol(Al s ^)] = flog" ^ - log" ±) < (log * V 

C 10g r b «T) Til C \ Tq Tq J TV. C \ Tq J 

as desired. 

ii) Let us perform the same changes of variables as in the previous part of the lemma. Then if 2 is the 
region in the Xi corresponding to W% and if $(^2) is the corresponding region in the yi from the condition 
specifying W2 we see that $(^2) this time is described by 

/ < log f o < Eti Vi < log i and Vi e [0, 00), 

{ < iog-^ < Ei=iw- 

Notice that the complement in A™ _r of the region $(f2 2 ) is described by 



aLh -$(n 3 ) 



< Ei=i 2/i < lo S 7^ and Vi € [0, 00 



, o<E-=i^<iog- 

Now introduce the unit Jacobean transformation 



so that finally we have 



3 

Zj = ^Vi, where j = 1.. 

i=l 



Vol($(ih)) = Vol(AV- n_) - Vol{complement) 

°g ro 



-^log"-^-- / / ... / / ... / <l:, ,<!-.„ 1. • • • il:.>d:.\ 



Zk-l J Zk 



k 1 / d 1 / „/n\ fe -J 



1 / R c/D\ n -^ 3 1 /. c/U 
'log log — log 



k— j j V 6 c ; 



As before from the bound on the Jacobean of $ we have 

VoZ(ft 2 ) < — VoJ($(fi 2 )) 
c 

and so we get the desired expression. □ 

2.1. Bounding the 7" term. On the set W\ we have a strong estimate on the sum in (|2.3p as follows. 
Notice that since < c/D < c we will have that the number of b's in the sums will be at least k, i.e. 
m > k, since all velocities corresponding to ro, Ti, ■ ■ ■ , r^ must have £ < 0. Let in fact b\, 62, ■ ■ ■ , 
be exactly the ones corresponding to ro, ri, . . . , r^ (or in other words to £0, £1, • • • , Cfe)- Observe that the 
requirement < c/D also ensures that for j < k + 1 

< 7T(c/£>) = 7T 

J - c-c/L> £>- 1' 

Furthermore notice that we always have at > tt, since if £ s corresponds to a% (and so r s > c, > 0) then 

7r|6, - c| r s 

Now recall that that I" is the same as the expression (|2.3p with replaced by Wi. Apply part ii) of 
Lemma 2.1. with B = -j^zi an d A = it to the J contribution (and bound &i = 7rro/|£ol by itself rather 
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than B which is certainly permitted in the proof of the Lemma) . Then after performing the integrations 
we have: 

I^<^(l + lM n n k (D-l)- k ^^ [ f[smaJ-Q-f + (t n )da n ...dai. 

If 1 I 4- 

Now we apply part i) of Lemma 2.2., recall that k = [fin] + 1 and bound < 1, /£" < M, to get 

^(l + l/nniD-lj/nr^ R'ro r 1 / R 
1 2tt r$ c \e \nl \ g r 

At this point note 2(1 + 1/V) < 2.7 and given the choice of f3 choose D so that 
(2.4) ((£> - l)/irf > 2(2.7) = 5.4 

which will give finally 

M i? 3 1/1 i? 
1 < 27rcro|f 1 |»l v 2 i0g r 
This clearly sums to an integrable function near £o = c (i.e. near ro =0), in fact we get 

M R 3b R 



E 

i=l 

For the pointwise estimate away from £ = c ? say f° r r o > c /2, see the remark at the end of this section 
which applies both to the terms Z™ and an d shows that there is no singularity near £q = 0. 
The convergence of the I\ term is also fast in the norm. If 7?.™ is the remainder term use 

(2.5) EtT^ 6 * 

i—n 

to see that 

MR' 1 r R |ej| / R 



Jr <R 2tt 2™n! 7 cr |Col \ r oJ 

?gration with respect to u — log ^ so dro — 

(2.6) / iej»< — - / u" e -»»d«< 

Jr„<fl c Z n - JO 



Switching to integration with respect to u = log ^ so <iro = —r^du with ro = i?e 11 we have: 

2Mi? 5 1 f°° „ s „ , 2Afi? 5 1 ,/2\" +1 4Mi? 5 



2™n! ' V 3 / 3™ +1 c 



2.2. Bounding the ZJ term. For this term we cannot obtain too strong a bound on the J contribution 
in (|2.3[) . but by Lemma 2.2. ii) we have a desirable bound on the measure of the integration region after a 
simple change of variables: namely that measure grows only as a fraction of n power of log ^ which will 
ensure that when we add the terms 1^ we will not get a nonintegrable singularity at ro = 0. 

As in the remark to Lemma 2.1. since J l ' m < J 1 ' 1 bound the J part of the integrand by -^b\. This 
remaining numerator term we can choose to be b\ where we arrange the bj so that b\ corresponds to fo so 
J < l/(27r)6i = 7rro/(27r|£i)|) and if we note |£i|/r„ < 1 we get 

(27r)™ 7rro i? 3 f i— l^ 1 1 2™ ro i? 3 /" ■". 

o n+2 o uT\ — T\sma i — f + {^ n )da n ...da 1 <—— r — T] sin ai da n ... da^ 

We now study the mass flux due to the remainder term TZ^ — Y)j^- n I\i i- e - we consider the remainder 
in L^(£o)- For this using part ii) of Lemma 2.2. we have 

J\Zo\<R q= ^ n - q plq - 4n Jo r o \€o\ c V c j V Dr J 

i=VM p R 



I? < 



2 n MR A 

7TC 

q—0,p—n — q 
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In the integral we make the substitution u = log so that dr a — —r du and Tq = j>e u and we get 



ro lo. 



— 

DrJ 



dr 



D 2 



q e~ 2u du 



„2 r0 



Ion 



D 2 



log- 



C 
1)2 



L> 2 



By noting that in A\ the integrand is maximized at the left end point and estimating the integral by the 
maximum of the integrand times the lenght of integration and then computing Ai explicitly we get 

M < (log £fj (^) 2 and M = [u^du = ^. 



With this we have 



l^|/ 2 "KoKo<™ 
i \<R ™ 



5 q=[Pn] 



' pig! \ c 



q— 0,p— n— <j ' 

The first sum is part of a binomial expansion and so 



"+1 nn»,n3 r 1=^ 



2 n MR A c 



E 



ttD 2 ^ n!29 

g— 0,p— n — g 



<?=[/8n] 

E 

g— 0,p— n — q 



l 



log 



DR 



n+l 



< 



log 



n+l 



The second sum is bounded by part of the remainder term in the Taylor explansion of ^-e 21og <= so as in 
(|2.5|) we have 



n!2 r ' 

g— 0,p— n — q 

with this we see 



v -^fiog^y<j- (2iog ^ ) 

^ «!5(n+i I s c/ 2™+! ((l-/3)n)! 



,2io g ^ _ 1 ( 21 °g- 



# 2 i? 2 



2 «+i ((l-p»! 



2™ M/? 5 2" 

ieS|jy«b)d& < i^iL^ 



log 



n+l 



ME 5 



2ttc ((1 -/?)«)! 



2 log 



(l-/3)n+l 



With the above bound on we interchange the integration and summation in the norm of the remainder 
7?-2 = X^n ^2 (which is permitted since the partial sums are bounded by an integrable function as pointed 
out in the singularity description further below). Thus we have 



L 



MR b log 



l£dl^2(£oRo < 



io\<R 



^U* / DRV ~ 1 



2 log 



DR 



(1-/3), 



C I 



In both sums above use that r! > r r e r . Thus in the second term use that ((1 — (3)i)\ > ((1 
/ g)j)(i-/3)i e -(i-/3)i so tna t this term is bounded by 



(2.7) 



2elog 



(1-/3)* 



< 



E 



-n+l 



e- 1 



provided that (2e log ^/((l - p») (1 ^ < 1/e. To ensure this recall that by (|2.4p we can choose D 
1 + 7r5.43 so log (RD/c) < 2 + 2//? + log(i?/c) and so (since i > n) it is sufficient to choose 

2e 1+ T^ 



n > 



1-/? 



„ 2 , R 
2 + -+log- 

p c 



For the first term similarly choosing n > 4e 2 (2 + 2/(3 + log — ) ensures 



(2. 



D -¥) 



< 



-n+l 



e- 1 
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Now noting that for n > 2 one has 3 „ 4 +1 < e ™ we combine the L]^ estimates on the remainders 1Z™ in 
(HU) and in ([23D and to get the desired result for K n = W{ + 11%, i.e. for 

n > 2 + - + log — max{4e^, — — } 



one has 



° c ) 1-0 

5/0 £?\ 



/ ^ oRo < 2 + - + log - e~ n+1 < 

J\£o\<R c \ [3 cj c 



ne n+1 . 



'\£o\<R 

If we now choose = 0.25 we get the statement of the theorem. We remark again that a more careful 
analysis may improve markedly the final constants. But for us it was most importnat that these constants 
are explicit and that we have established the convergence mechanism. This completes the estimates of the 
mass flux of the remainder terms in our theorem. 

Next we show the singularity description. Using part ii) of Lemma 2.2. where we bound the terms in 
the binomial expansion by the highest power in which they appear we have that 
->3 r on n j i?3 



2 n r R 6 f A . 2 n M R 3 r r 2" 





c/D 
log 






ro 



0n 



With this we get that the total contribution of the terms is 

log 



M R 3 ^ 1 ,„ A RD 



=i 



c/D pn 



ro 



(2.9) h < ^-^log^/^llog^/Dro)!^^^'^! 1 -^!^, 

where Ci(D,c,R) = 41og^. This is an integrable function near r = since as < 1 for all S > we 
have 

i.e. the term with the exponent grows slower than any negative power of ro. With this in mind, only for 
the purpose of the singularity description, we pick on the right hand side of (|2 .4(1 instead of the choice 
there (so that in the J™ term the exponent of r is —(1 + S) instead of —1.5). So together with (|2.9[) we 
complete the singularity description in part ii) of our main theorem: 

MR 3 1 

/-(Co) =h + h< C s (c,R) 1+s , 

c lsol r o 

provided ro < c/2 and where the constant Cs(R,c) is explicitly derived from the expression (|2.9p . When 
ro > c/2 we note that the remark at the end of this section further shows that there is in fact no singularity 
near £5- = 0. 

With this the proof of our theorem is completed. 

Remark 2.4. For the pointwise estimate of / - (£o) we consider the case ro > c/2 separately because 
in the above we have b\ for a remaining term in the sum and b\ = is unbounded for |£q| near 

(which can happen if ro > c/2). This point is very easy to correct however since now we are not looking to 
optimize the singularity at ro = 0. In fact we will bound the sum part of the integrand exactly as in the 
previous section, except this time we choose the term remaining in the denominator to be a\ rather than 
b\ where we agree that a; is in fact the term corresponding to £„, i.e. ai = j c ^ . Notice that then this 

If 1 1 

term cancels out with the factor in front of /(£„) in the formula. Thus without splitting the set W 
this time and just using an argument as the one in part i) of Lemma 2.2. (where the lower bound on the 
sum of the y\ is now just so we simply bound the Jacobean of <& by 1) we get for ro > c/2 the bound: 

(27t)" R 3 f A • » r , , 1 i? 3 2" / R Y , , 

fn < y**^ y^nsma^...^ < ^-j M 
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The sum of these terms is relevant for r > c/2 and has no singularity at £q = and vanishes at r = R 
which is the boundary of the domain of influence. ■ 

Remark 2.5. One would hope to extend the above proof to the case when the background is a Gaussian 
distribution as is very natural for problems in kinetic theory. It is expected that the same basic convergence 
mechanism will exist in this case but technical difficulties will arise due to the fact that K now will not 
have the strict shrinking property. In other words due to a diffusive component (in magnitude related to 
the temperature) the sequence r, from the proof does not have to be increasing - however note that this 
sequence will be close to increasing with a very high probability which may lead to a very similar situation. 

■ 

3. Appendix - Carleman's representation for the gain term 
Consider the gain term 



Q + (f, 9)(0= / / f(£-((S-t.)-n)n)9(£. + ((t-t.)-n)n)\(£-£.)-n\d£.dn 
Jr 3 Js+ 

which according to the geometric intuition explained in the introduction we would like to write as 

Q + (/,.9)(£o)= f g&)W g &,So) f f{QW f (t.,So,Od£dZ t 
Jr 3 Jl ( , (0 

where W g and Wf are weight functions to be determined. We will use the symmetries of the '-transformation 
so we observe 

/ <KQQ + {f,9)(tw= f [ /(«*)/ ^em-^)-n\dnd^dc 

JR3 JR3 JR3 Js+ 

Now the velocities form the vertices of a rectangle and with £ and fixed this rectangle is 

inscribed in a fixed sphere in different ways according to the value of n (n is parallel to — 

If da is the surface measure on that sphere one easily finds that 

\(£ - ■ n\dn = JjTZJ^- 

Therefore 

I mQ + {f,9)im= I [ / <t>(o jr }- ri d*d£ t d£. 

JR3 JR3 JR3 JS(t,£,) It - S*l 

and we will take <fi e ((,) — ^ttI) 3 ' 2 e~^°~^ 2 /( 2e ) so that we can compute 

Q + (.f,g)(to) = lim f M0Q + tf,g)(0dt 

After denoting z to be the midpoint between £ and (which is the center of the sphere for the da integral) 
a straightforward calculation using the above representation shows 

MOQ + (f,g)(Od{ = 

(2ire) L / z \z - £ | L J 
Notice this symmetric expression shows that Q + (f, g) = Q + (g, /), but of course Q(f, g) ^ Q(g, /) due to 
the nonsymmetric loss term. It is obvious that the second term in the brackets above does not contribute 
in the limit e — > 0. Then a careful but straightforward calculation shows that the first term in the limit is 
equivalent to j^z^S(h) where h is the distance from the plane i^,^ and we indeed arrive at the desired 
formula 



Q + (f,g)(to) = / gi&rr^-n I /Km.- 

JR3 It* — ?0| JL itiQ 
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